Reducing decoherence in optical and spin transitions in rare-earth-ion doped materials 
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In many important situations the dominant dephasing mechanism in cryogenic rare-earth-ion 
doped systems is due to magnetic field fluctuations from spins in the host crystal. Operating at a 
magnetic field where a transition has a zero first-order- Zeeman (ZEFOZ) shift can greatly reduce 
this dephasing. Here we identify the location of transitions with zero first-order Zeeman shift 
for optical transitions in Pr 3+ :YAG and for spin transitions in Er 3+ :Y2Si05. The long coherence 
times that ZEFOZ can enable would make Pr 3+ :YAG a strong candidate for achieving the strong 
coupling regime of cavity QED, and would be an important step forward in creating long-lived 
telecommunications wavelength quantum memories in Er 3+ :Y2SiOs. This work relies mostly on 
published spin Hamiltonian parameters but Raman heterodyne spectroscopy was performed on 
Pr 3+ :YAG to measure the parameters for the excited state. 

PACS numbers: 03.67.Pp, 76.30.Kg, 76.70.Hb, 42.50.Md 



I. INTRODUCTION 

In a rare-earth-ion-doped system there is a contribu- 
tion to the homogeneous linewidth from the interaction 
between the dopants and the atoms that make up the 
host crystal Fluctuations in the electronic and nu- 
clear spin of the host atoms cause the dopants to expe- 
rience small random magnetic fields that can alter the 
transition frequency via the Zeeman effect. This has the 
effect of broadening the transition. In systems with nu- 
clear spin greater than one applying a magnetic field to 
the sample produces anticrossings in the energy of the 
4^ . levels, and around these anticrossings it is possible to 
find transitions that have a zero first-order Zeeman shift 
[2|45|. This means that small fluctuations in magnetic 
field will no longer change the transition frequency, thus 
decreasing the broadening associated with host-dopant 
interactions. 

The technique of using a ZEFOZ transition to negate 
the ion-spin interaction was first used by Fraval et al. 
in an attempt to lengthen the hyperfine coherence time 
in Pr 3+ :Y2Si05, as it was deemed too short for prac- 
tical quantum computing applications. They initially 
demonstrated an increase in phase memory time from 
550 /is to 82 ms [2]. By improving the stability of the 
setup this time was increased to 860 ms Q . The correla- 
tion time for the remaining frequency perturbations was 
long and the same work used dynamic decoupling tech- 
niques to increase echo decay time coherence beyond 10 s. 
The ZEFOZ technique has also been used to measure ex- 
tended coherence times in Pr 3+ :La2(W04)3 [Bj], and the 
locations of the ZEFOZ transitions for Eu 3+ :Y2Si05 and 
EUCI36H2O have been determined Q. 
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In this work we provide an extension to the ZEFOZ 
technique by showing that it can be applied to opti- 
cal transitions and the hyperfine transitions of Kramers' 
ions. 

Reducing the decoherence of an optical transition has 
previously been investigated in Pr 3+ :LaF3 by Macfarlane 
et al. [6]. By driving the nuclear spins of the fluorine 
atoms with a radio- frequency (RF) field they showed that 
it was possible to decouple these spins from the Pr 3+ ions. 
By choosing the RF field strength such that the angle of 
the effective field experienced by the fluorine atoms is 
the 'magic angle', the coherence time was increased by 
an order of magnitude. The ZEFOZ technique should be 
more effective at reducing the decoherence because it will 
reduce the interaction of the dopants with all of the host 
atoms as opposed to just one type of atom. 

As an example, the ZEFOZ transitions are found for 
the optical transition at 609.76 nm in Pr 3+ :YAG. Based 
on the measured values of the population lifetime and 
coherence time in Pr 3+ :YAG, it should be possible to 
increase the optical coherence time by an order of mag- 
nitude. If the coherence time can be increased, it will 
make Pr 3+ :YAG an excellent candidate for achieving the 
strong-coupling regime in cavity QED |7]. To find a ZE- 
FOZ transition in Pr 3+ :YAG, knowledge of the entire 
spin Hamiltonian, for both the ground and excited states, 
is required. Here the necessary spin Hamiltonian parame- 
ters are measured using Raman heterodyne spectroscopy. 

To date all applications of the ZEFOZ technique have 
focused on non-Kramers' ions; however, it would also be 
useful to determine whether it is possible to find ZE- 
FOZ transitions for Kramers' ions. As an example, the 
location of hyperfine ZEFOZ transitions have been cal- 
culated for Er 3+ :Y 2 Si05. Er 3+ :Y 2 Si0 5 is an interesting 
system to study because it has the longest optical coher- 
ence time of any rare-earth- ion system [§| . If a hyperfine 
transition could be found that also had a long coher- 
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ence time, it would make Er 3+ :Y2Si05 an ideal system 
for creating a quantum memory at telecommunication 
wavelengths. Here the location of ZEFOZ transitions in 
g r 3+ : Y 2 Si05 are calculated using the spin Hamiltonian 
parameters measured by Guillot-Noel et al. @. 



II. HYPERFINE INTERACTION IN 
RARE-EARTH-IONS 

The 4f energy levels of a rare-earth-ion are completely 
described by the following Hamiltonian [lQj |: 

H = H FI + Hcf + Hhf + Hq + H eZ + H n z, (1) 

where the first two terms represent the free-ion and 
crystal-field Hamiltonians and determine the electronic 
energy level structure. The other four terms are much 
smaller than the first two and describe the hyperfine 
structure of the energy levels. These terms are the hyper- 
fine, quadrupole, electronic Zeeman and nuclear Zeeman 
Hamiltonians respectively. As the free-ion energy levels 
are well known, and only the lowest crystal-field level is 
populated at cryogenic temperatures, only the final four 
terms of Eq. (pQ) need to be considered when calculating 
the energy of hyperfine levels. 

In this work two different materials are considered; 
Pr 3+ :YAG which is a non-Kramers' ion (has an even 
number of f-electrons), and Er 3+ :Y2Si05 which is a 
Kramers' ion (odd number of f-electrons). Different 
Hamiltonians are used to describe the two systems. 

For a non-Kramers' ion in a site of low symmetry the 
electronic angular momentum is quenched, lifting the de- 
generacy of the electronic states. Applying second-order 
perturbation theory to the system causes second-order 
hyperfine and second-order electronic Zeeman terms to 
arise that are of a similar magnitude to the nuclear 
Zeeman and quadrupole interactions. Therefore, non- 
Kramers' ions can be described by the effective Hamilto- 
nian derived by Teplov [TTj |: 



# eff = gjfi B B A B + B (2Ajgjfi B A + hj N ) • J 
+ /-(A 2 J A + T Q )-/, 



(2) 



where gj is the Lande g value, hb is the Bohr magne- 
ton, 7at is the nuclear gyromagnetic ratio, and Aj is the 
magnetic hyperfine constant. B is the applied magnetic 
field, J is a vector of the nuclear angular momentum op- 
erators, and A is a tensor that describes the second-order 
interactions in the system. 

The term g 2 jji B B K B describes the quadratic elec- 
tronic Zeeman interaction [l0|. B ■ (2Ajgjjj jB A-\-h / yN) ■ I 
is the enhanced nuclear Zeeman interaction, where the 
nuclear Zeeman interaction is enhanced due to cou- 
pling to the second-order hyperfine interaction fl2l [l3| . 
/ • (A 2 jA + Tq) • / is the combined quadrupole inter- 
action, made up of the second-order hyperfine interac- 
tion and the pure-quadrupole interaction. Because the 
second-order hyperfine interaction can be written in the 



same form as the pure-quadrupole interaction, it is often 
called the pseudo-quadrupole interaction fl4l |. 

For the work done on Pr 3+ :YAG the dopants occupy 
sites of sufficiently high symmetry that the axes of the 
interaction tensors coincide with the site axes. This en- 
ables the Hamiltonian to be diagonalized, i.e. 

H e fl = BZB + BMI + IQI, (3) 

where M is the enhanced nuclear Zeeman tensor: 
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Q is the tensor of combined quadrupole constants: 
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and Z is the second-order Zeeman tensor: 



Z = 











z z 



When considering ions with an odd number of electrons 
(such as Er 3+ ), Kramers' theorem states that all levels 
will have an electronic degeneracy that can only be lifted 
by applying a magnetic field to the sample |10(. This 
results in a large contribution to the Hamiltonian from 
the first-order electronic Zeeman and magnetic hyperfine 
interactions. At cryogenic temperatures only the lowest 
crystal-field level is populated; therefore, the system can 
be approximated as having spin S = \ Q . The enhanced 
nuclear Zeeman and combined quadrupole interactions 
are as for non-Kramers' ions. Thus, the energy levels of 
Kramers' ions can be described by the following effective 
Hamiltonian [Ioj |: 

tfeff = VbB g S + I A S + I (A 2 jA + T Q ) • J 

- B ■ (2Ajgjii B A + Hj N ) ■ J, (4) 

where \±bB • g ■ S is the electronic Zeeman term, and 
I - A - S is the first-order hyperfine interaction. 

To simplify the description for Er 3+ :Y2Si05 we will 
rewrite the above Hamiltonian as [9(: 

tfeff = n B B - g. S + I- A- S + I- Q.I - fi n g n B • I, (5) 

where ji n is the nuclear magneton and g n is the nuclear 
g-factor = -0.1618 for 167 Er . g, A, and S are the g- 
factor matrix, matrix of hyperfine parameters and the 
electronic spin matrix respectively. 



III. PR 3+ :YAG 

Y3AI5O12 (YAG) has a cubic crystal structure. The 
praseodymium ions substitute for yttrium ions that are 
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located at do decahedral c-sites, and have local D2 sym- 
metry flEl [~L6| . This means there are six different ori- 
entations that the praseodymium ions, which substitute 
for yttrium, can have in the crystal. This is explained in 
more detail in the context of Tm 3+ :YAG in [16]. The six 
sites are magnetically inequivalent, meaning the ions' hy- 
perfine levels will experience a site-dependent frequency 
shift under the application of an external magnetic field. 
The orientation of the local a>, y-, z-axes of site 1 are 
defined as being along the [1,1,0], [1,-1,0] and [0,0,1] di- 
rections of the crystal. The orientations of the other sites 
are found by performing an appropriate rotation on the 
site 1 axes. 

The transition of interest is between the X D2 excited 
state and 3 H4 ground state in Pr 3+ :YAG, occurring 
at 609.76 nm (in vacuum) (l7| . The inhomogeneous 
linewidth of this transition has been measured as 1 - 
2 cm -1 [TH, [13, EE] and has an oscillator strength of 
/ = 1.3 x 10 -6 pjj. Shelby et al measured an opti- 
cal coherence time of T 2 ( pt) = 20 [is [l8[ compared to an 
excited state lifetime of Ti( opt ^ = 230 /is pjj. They put 
the discrepancy down to magnetic dephasing caused by 
spins in the host. The coherence can be expected to be 
shorter for Pr 3+ :YAG than for example Pr 3+ :Y 2 Si05 be- 
cause the nuclear magnetic moment of aluminium (3.64 
(2p|). the dominant spin in YAG is much larger than for 
yttrium (-0.137 [20|), the dominant spin in Y2Si05. 

If this excess dephasing is purely magnetic, it should 
be easily reduced to a level where the coherence time is 
lifetime limited, requiring only a twenty fold increase in 
coherence time to reach a T2 = 2T\. This can be com- 
pared to the thousand fold increases demonstrated for 
hyperfine transitions. Often in rare earth doped systems 
the broadening due to phonon scattering is negligible at 
cryogenic temperatures [loj]; however, there is a possi- 
bility that this could be a more significant issue with 
Pr 3+ :YAG because of a crystal field level not far from 
the ground state [15j. 

The hyperfine structure of Pr 3+ :YAG was first studied 
by Shelby et al. where spectral hole burning was used 
on a 0.15% Pr 3+ :YAG crystal to measure the hyperfine 
splittings of the ground and excited states [l8|. Pr 3+ has 
zero electronic spin and a nuclear spin of 5/2, resulting 
in three degenerate pairs of hyperfine levels for both the 
ground and excited states. Wang later confirmed these 
hyperfine splittings using optically detected nuclear mag- 
netic resonance [15] , and they also agree with what we 
have measured using Raman heterodyne spectroscopy. 

An in-depth investigation into the hyperfine structure 
of Pr 3+ :YAG was performed by Wang [15]. Based on the 
spin Hamiltonian of the system and the measured hyper- 
fine splittings, he determined the values of the effective 
quadrupole tensor for both the ground and excited states. 
Measurements were also performed using Raman hetero- 
dyne spectroscopy on the 33.33 MHz ground- state hyper- 
fine transition. By observing how the transition energies 
varied with magnetic field, the values of the enhanced 
nuclear Zeeman tensor were calculated. The known spin 



TABLE I: Spin Hamiltonian parameters measured by Wang 



[15] 



Ground State ( 3 H 4 ) Excited State ( X D 2 



5x (MHz) 
S 2 (MHz) 
D (MHz) 
E (MHz) 
lx (kHzG- 1 ) 
ly (kHzG- 1 ) 
lz (kHzG" 1 ) 



33.33 
41.86 
11.36 
2.776 
13 
1.6 
30.3 



6.41 
8.32 
2.24 
0.520 
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FIG. 1: (Color online) Example Raman heterodyne spectrum 
centred on the 33.33 MHz ground-state transition. A beat is 
detected when the RF field applied to the sample is resonant 
with a hyperfine transition. 



Hamiltonian parameters are summarized in Table |TJ 

To find an optical ZEFOZ transition in Pr 3+ :YAG, 
knowledge of the entire spin Hamiltonian for both the 
ground and excited states is required. Therefore, the 
excited-state enhanced nuclear Zeeman tensor needs to 
be measured. This can be done using Raman heterodyne 
spectroscopy [H|,[22j]. 



A. Measuring Spin Hamiltonian Parameters 

Raman heterodyne spectra were obtained using an ex- 
perimental setup similar to that described in reference 
|23j . An upgraded set of Helmholtz coils was used that 
enabled DC magnetic fields of up to 200 G to be gen- 
erated. When recording the spectra the sample was im- 
mersed in super-fluid liquid helium to produce tempera- 
tures less than 2.2 K as this gave the strongest Raman 
heterodyne signal (most likely due to the sample being 
colder) . Figure [1] is an example of the recorded spectra. 



The crystal used in the experiments was a 0.5% 
Pr 3+ :YAG crystal with an initially unknown orienta- 
tion, grown at the Australian National University. The 
crystal was placed in the cryostat so that the crystal 
faces are perpendicular to the axes of the Helmholtz 
coils. Using Raman heterodyne spectroscopy on the 
33.33 MHz ground-state transition (the same transition 
used by Wang to measure the ground-state parameters), 
the crystal orientation was determined to be that the 
Y-, Z-axes are approximately aligned along the [1,-1,0], 
[1,1,1] and [1,1,-2] directions. 

Raman heterodyne spectroscopy was then performed 
on the 6.41 MHz transition so that the components of 
the excited-state enhanced nuclear Zeeman tensor can be 
determined. The recorded spectra are shown in Fig. [2j 
Each plot corresponds to a series of spectra (each one 
similar to Fig. [1]) laid side by side, each with a different 
DC magnetic field applied to the sample. The magnetic 
field coils perform a rotation of 180° around the X-Z 
plane, 180° around the X-Y plane, and 180° about the 
Y-Z plane (X, Y, Z directions are defined in terms of the 
three sets of Helmholtz coils). Due to the weakness of this 
transition, only four of the twenty-four possible lines were 
visible. Each measured spectrum was inspected manually 
and the frequencies of the peaks corresponding to tran- 
sitions were selected. These frequencies are represented 
by the black crosses in Fig. [2j 

To calculate the excited state nuclear Zeeman parame- 
ters a computer program was run that minimized the mis- 
fit between the measured data and the theoretical values 
(calculated from the spin Hamiltonian). There were two 
extra factors that needed to be taken into consideration 
when performing the fit. Firstly, an extra rotation was 
applied to the spin Hamiltonian to allow for imperfect 
alignment between the sample and the coil axes. This 
extra rotation is defined as: 



1 

cosa sina 
— s'ma cosa 



cos/3 —sin/3 

1 
sin/3 cos/3 



COS7 sin7 
— SU17 COS7 
1 



The second consideration was that the principle axes of 
the combined quadrupole tensor are not necessarily the 
same for the ground and excited states. However, due 
to the local D2 symmetry of the praseodymium ions, 
the assumption was made that the axes of the com- 
bined quadrupole tensor should lie along the three Co- 
axes ([1,1,0], [1,-1,0], and [0, 0, 1] for site 1). The best fit 
was found for: 
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FIG. 2: Raman heterodyne spectra obtained for the excited- 
state hyperfine transition in Pr 3+ :YAG centred on 6.41 MHz. 
The horizontal axis corresponds to the orientation of the ap- 
plied DC magnetic field (Bo = 68 =b 2 G) , the vertical axis 
the frequency of the applied RF, and the shades of gray rep- 
resents the intensity of the beat signal. Black crosses are the 
peaks selected from each individual spectrum. Solid lines are 
the best fit to the crosses. The spectra correspond to the 
magnetic field being rotated about the; (a) X-Y plane, (b) 
Y-Z plane, and (c) X-Z plane. 
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and is represented by the solid lines in Fig. [2j The pa- 
rameters of this fit are given in Table [TTJ 

When fitting to the Raman heterodyne data the 
second-order Zeeman term (B • Z • B) was excluded from 
the calculations. While this term adds a frequency shift 



dependent on the magnitude of the applied magnetic 
field, the shift is constant for all hyperfine levels of a 
particular state. Therefore, the frequency of the hyper- 
fine transitions is independent of this term. However, it is 
necessary to include the second-order Zeeman term when 
calculating the relative frequencies of optical transitions. 
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TABLE II: Parameters of the best fit to the data obtained for 
the 6.41 MHz excited-state transition of Pr 3+ :YAG. 



TABLE III: Components of the second-order Zeeman tensor 
in Pr 3+ :YAG. 



Excited- State Parameter Value 



lx (kHzG -1 ) 
lv (kHzG" 1 ) 
lz (kHzG -1 ) 
a 

P 

7 



1.39 
1.68 
1.62 
-2.7° 
0.4° 
3.8° 



Ground State ( 3 H 4 ) Excited State ( X D 2 



9J 

(J\\N\\J) 
Aj (MHz) [24] 
Z x (kHzG -2 ) 
Z y (kHzG -2 ) 
Z z (kHzG -2 ) 



0.8 

1.316 

1093 
9.54 x 10" 4 
2.41 x 10" 5 
2.37 x 10" 3 



1.0 
1.000 
831 
1.15 x 10" 5 
5.03 x 10" 5 
4.37 x 10" 5 



B. Finding the ZEFOZ transitions 

The motivation behind measuring the spin Hamilto- 
nian parameters for Pr 3+ :YAG is to find an optical tran- 
sition where there is zero first-order Zeeman shift, result- 
ing in an extended coherence time. Before calculating 
what magnetic field to apply to the sample, it is nec- 
essary to calculate the components of the second-order 
Zeeman tensor, Z = g 2 jH 2 B K. The second-order hyper- 
fine tensor is calculated from second-order perturbation 
theory [Hj: 



2J+1 



n=l 



(0| J a \n) (n\ Jp |0) 
E n — Eq 



(6) 



where |0) is the level of interest, \n) refers to the other 
crystal field levels, and E is the energy of the level. Here 
A is calculated by subtracting the nuclear gyromagnetic 
ratio (7at = 1.304 kHzG -1 for praseodymium) from its 
enhanced value (j a = j x ,Jy,Jz), that is: 



A = (la ~ In) h 
2g J /j jB A J 



where the magnetic hyperfine constant is given by: 

Aj = 2^ B7N h(r- 3 }(J\\N\\J). (8) 

(r -3 ) is the expectation value of the inverse-cube 
electron- nuclear distance ((r -3 ) = 5.0 a.u. for Pr 3+ j2~" 
and: 



(J\\N\\J) 



1 



J(J+1) 
- 3(L ■ J)(L ■ S)]}, 



[(L-J) + £[L(L+1)(S-J) 



(9) 



21+1-45 



-. In this work the values of 



where £ = S (2e-i)(2£+3)(2L-i)' 
the magnetic hyperfine constant given in reference [24[ 
are used. These were calculated from magnetic resonance 
measurements on salts, but should also be applicable to 
crystalline hosts. 

From Eqs. ([7j)-([9]) the parameters of the second-order 
Zeeman tensor are calculated. These are summarized in 
Table HIl 



When calculating A it was necessary to make certain 
assumptions about the sign of the enhanced nuclear gy- 
romagnetic ratios. Raman heterodyne spectroscopy only 
determines the magnitude of the enhanced nuclear gy- 
romagnetic ratios, as inverting the principle values of 
the enhanced nuclear Zeeman tensor leaves the spectra 
unchanged. This effectively means A aa = k\j a — jn\ 
(k = h/2gj/j J BAj), which has two solutions. For the ex- 
cited state the enhanced nuclear gyromagnetic ratios are 
approximately the same size as 7at; therefore, either A aa 
is small, or A aa ~ 2kjN. It is assumed that the former is 
more likely; therefore, the sign of j a was chosen such that 
1 7a — 7aH < 7a- The situation is different for the ground 
state in that ^ x and j z are significantly larger than jjy. 
In this case it is difficult to make an assumption about 
the sign of 7^; however, either choice will result in similar 
values of A acn that is, A aa ~ k\j a \. Here the sign of j a 
was again chosen such that \j a — jn\ < j a . It is thought 
that the resulting landscape of ZEFOZ transitions will be 
similar for both choices of 7^, because for the magnetic 
fields at which ZEFOZ transitions exist (10 2 - 10 4 G) the 
magnitude of B • Z ■ B is smaller than the magnitude of 
both B • M • / and I Q I. Therefore, even with the 
wrong choice of j a , it should be possible to find the ZE- 
FOZ transitions experimentally by tuning the magnetic 
field. 

If the ZEFOZ transitions cannot be found experimen- 
tally there are two options. The first is to calculate the 
ZEFOZ transitions for the alternative values of A aa . The 
second, better option would be to measure the second- 
order Zeeman tensor, and recalculate the ZEFOZ transi- 
tions using this information. The second-order Zeeman 
tensor could be measured by burning a spectral hole, then 
measuring the frequency shift of the hole when a mag- 
netic field is applied. 

Now that the entire spin Hamiltonian is known, for 
both the ground and excited states, the optical ZEFOZ 
transitions can be found. This is done by adopting a 
method similar to Longdell et al. [H . By searching over 
a three-dimensional grid of magnetic field values, tran- 
sitions are found with a gradient vector equal to zero. 
A three-dimensional Newton- Raphson method is used to 
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FIG. 3: (Color online) Magnitude of the applied magnetic 
field versus the maximum curvature of the ZEFOZ transitions 
found for Pr 3+ :YAG. 



hone in on the zero points. The iteration: 

Bnew — B 



V 

2C' 



(10) 



is used, where v and C are the Zeeman gradient and 
curvature tensors respectively: 



^(B) = ^(B)- U ,(B)], 



CfHB) 



d 2 



dBidB, 



■[w p (B)-w,(B)]. 



(11) 



(12) 



These are calculated using first- and second-order per- 
turbation theory: 

^- Wp (B) = (^(B)|C y -|^(B)), (13) 



d 2 



dBidB 



■Wp(B) 



/> p (B)|q|^(B))(^(B)|^|^(B)) 
u p (B) - u q (B) 



(14) 

where ZijBj + BiZij + M^/y. The ZEFOZ transi- 

tions found for Pr 3+ :YAG are shown in Fig. [3l A total 
of 122 ZEFOZ transitions were found with curvatures 
ranging from 10 1 to 10 5 HzG -2 . The magnitude of the 
required magnetic fields ranges from 50 to 10000 G. The 
transition with the lowest curvature should result in the 
largest increase in coherence time. This is because hav- 
ing a lower curvature means the transition frequency is 
more resilient to small fluctuations in magnetic field. 

In previous ZEFOZ experiments, operating on a turn- 
ing point with curvature of 60 HzG -2 resulted in a 1700- 
fold increase in the coherence time of Pr 3+ :Y2Si05 
likewise, in Pr 3+ :La2(W04)3 a turning point with curva- 
ture of 120 HzG -2 resulted in a 630-fold increase [5|. In 



Pr 3+ :YAG only a 23-fold increase is required to achieve 
a population lifetime limited coherence time (from 20 /is 
to 460 /is), and turning points exist with curvatures as 
low as 10.6 HzG -2 . Therefore, by operating on a ZEFOZ 
transition we expect to see a significant increase in the 
optical coherence time. 



IV. ER 3+ :Y 2 SI0 5 

There are two main factors that make Er 3+ :Y2Si05 
an interesting material system for studying with respect 
to rare-earth-ion doped quantum computing. Firstly, 
the 4 Ii5/2 - 4 Ii3/2 transition occurs at 1536 nm (site 
1) and 1538 nm (site 2). This is important because it 
falls within the telecommunications C band. For exam- 
ple, were a working quantum repeater to be realized from 
g r 3+ : Y 2 gi05, existing telecommunications infrastructure 
could be used, greatly reducing the cost of implementing 
a quantum network. Secondly, the 4 Ii3/2 excited state 
of Er 3+ :Y2Si05 has a long optical population lifetime 
(Tu npt) = 11.4 ms for site 1, Ti{p P t) — 9-2 ms for site 2 
|25l|), and the longest optical coherence time (T 2 f op t) = 
4.38 ms jH) of any solid-state system. The aim of this 
work is to find a ZEFOZ transition in Er 3+ :Y2SiOs, with 
the intention of extending the coherence time of a hyper- 
fine transition. This should allow long-term storage of 
quantum states for quantum computing applications. 

There have been several investigations into the lifetime 
of the ground state levels in Er 3+ :Y 2 Si0 5 [H,^. Using 
spectral hole burning Hastings-Simon et al. measured 
the ground state lifetime between electronic Zeeman lev- 
els to be T^hypj = 130 ms; however, they also observed 
that a subset of the ions had lifetimes up to 60 s (26| . 
Baldit et al. measured the hyperfine lifetime of isotopi- 
cally pure 167 Er to be Tuh yp \ = 97 ms |27[. Based on 
the inequality T2 < 2T\ it should be possible to attain 
hyperfine coherence times up to T 2 (h yp ) = 194 ms. 

Erbium has six stable isotopes of which only 167 Er has 
a non-zero value of nuclear spin. This isotope is of par- 
ticular interest as the non-zero nuclear spin results in 
energy levels that exhibit hyperfine structure. 16T Er has 
a natural abundance of 22.95% [28| and a nuclear spin of 
/ = I . For the J = ^ ground and J = ^ excited states 
the crystal field splits into eight and seven Kramers' dou- 
blets respectively. As stated earlier the system is approx- 
imated as having spin S = \ ; thus each crystal- field level 
will split into a total of sixteen hyperfine levels (two elec- 
tronic spin states (ms = |, ^7-) x eight nuclear spin 



states (mi = ^ - |)). 

Using electron paramagnetic resonance spectroscopy, 
g, A, and Q have been measured by Guillot-Noel et al. 
for the 4 Ii5/2 ground state of sites 1 and 2 [9j. Sun et al. 
used optical Zeeman spectroscopy to measure g for the 
ground and excited states of both sites [29j. These two 
separate measurements of g for the ground state agree ex- 
tremely well, except for along the principle x-axis of site 
1. This discrepancy is because the magnitude of g in the 
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x direction is significantly smaller than in the other di- 
rections; therefore, the measurements are less sensitive to 
variation along this axis. For the calculations performed 
in this work, the parameters measured by Guillot-Noel 
et at. [9] are used as they constitute a complete mea- 
surement of the ground-state spin Hamiltonian. For site 
1: 



2.92 -3.08 -3.68 
-3.08 8.19 5.96 
-3.68 5.96 5.52 



Ax 



69.35 -580.73 -248.83 
-580.73 696.30 682.49 
-248.83 682.49 495.54 



MHz, 



Qi 



21.40 -8.18 
-8.18 3.79 
-15.27 0.60 



-15.27 
0.60 
-25.20 



MHz, 
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and for site 2: 



g 2 



14.75 -2.02 
-2.02 1.89 
2.62 -0.93 



2.62 
-0.93 
0.05 



-1521.40 178.11 -141.76 
178.11 172.09 212.54 
-141.76 212.54 199.01 



MHz, 



FIG. 4: (Color online) Energy of the site 2 hyperfine levels 
as a magnetic field is applied along the Di-axis of the crystal 
for; (a) low magnetic fields, and (b) high magnetic fields. To 
gain an understanding of the relative frequency between levels 
at high fields, in (b) the frequency of the lowest eight energy 
levels relative to the 4th lowest energy level are plotted. For 
clarity an extra offset has been added to each level such that 
the relative frequency is zero at B = 1 kG. (Note that both 
the x and y axes for the graphs are different.) 



Q 2 = 



-3.50 -19.84 
-19.84 50.40 
24.22 6.73 



24.22 
6.73 
-46.90 



MHz. 



The coordinates are defined in terms of the axes of the 
Y2Si05 crystal, that is, x = Di-axis, y = 1^2-axis, z = b- 
axis. 

Knowledge of the entire spin Hamiltonian allows the 
frequency of the hyperfine levels to be calculated for any 
applied magnetic field. As an example, Fig. H]^a) shows 
how the ground-state hyperfine splittings change as a 
magnetic field is applied along the Di-axis of the crystal. 
Notice that the energy levels exhibit a large number of 
anticrossings; therefore, ZEFOZ transitions should exist. 
At large applied fields the electronic Zeeman interaction 
dominates and the system displays two distinct levels cor- 
responding to the two electronic states, each of which is 
split into eight hyperfine levels. 

To gain an understanding of why ZEFOZ transitions 
exist even for large magnetic fields the relative frequency 
of the wis — — 1/2 hyperfine levels is plotted in Fig.|4^b). 
The relative fluctuations in transition frequency are small 
(on the order of a few MHz) which will result in ZEFOZ 
transitions with a low curvature. 



A. Finding the ZEFOZ transitions 

Similar to the earlier work on Pr 3+ :YAG, a three- 
dimensional Newton- Raphson method is used to find the 
hyperfine ZEFOZ transitions in Er 3+ :Y2Si05. The gradi- 
ent and curvature are as defined in Eqs. (pT]) - ([Tl)) except: 



Cij 



V>n9nh- 



(15) 



Figure [5] is a plot of the absolute value of the applied 
magnetic field versus the maximum value of the curvature 
for the identified ZEFOZ transitions. For site 1 a total 
of 102 ZEFOZ transitions were found, with curvatures 
ranging from 10 1 to 10 8 HzG -2 and requiring magnetic 
fields ranging from 0.1 to 100 kG. 635 ZEFOZ transitions 
were found for site 2, with curvatures ranging from 10 -2 
to 10 8 HzG- 2 . The required magnetic fields range from 
0.01 to 100 kG. 

The maximum curvature of the best site 2 ZEFOZ 
transitions, were three orders of magnitude smaller than 
for those of site 1. Therefore, when experimentally trying 
to measure extended coherence times it is sensible to fo- 
cus on site 2, as the transition with the lowest curvature 
should give the largest increase in coherence time. 
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FIG. 5: (Color online) Magnitude of the applied magnetic 
field versus the maximum curvature of the ZEFOZ transitions 
found for; (a) site 1, and (b) site 2 of Er 3+ :Y 2 Si0 5 . Note 
the two shaded regions in (b) which correspond to the two 
different regimes as described in the text. 



B. Properties of the ZEFOZ transitions 

In Fig. [5](b) there are two different shaded regions. 
The shaded green (right) region corresponds to ZE- 
FOZ transitions in the 'strongly mixed' regime, whilst 
the shaded red (left) region corresponds to the 'weakly 
mixed' regime. The strongly mixed regime is where the 
electronic Zeeman splitting (the term /i^B-g-S in Eq. [5J 
is either small or comparable to the hyperfine splitting 
(I • A • S). The weakly mixed regime is where the elec- 
tronic Zeeman term dominates and the electron spin is 
close to a good quantum number. Figure 4 shows the 
energy levels for site 2 as the field is increased along the 
crystal .Di-axis. For low fields (< 300 G for this field 
direction) there are many avoided crossings due to the 
interplay between the hyperfine splitting and the elec- 
tronic Zeeman term. It is in this area where the strongly 
mixed ZEFOZ transitions occur. At high fields where the 
electronic Zeeman term is larger than the hyperfine split- 
ting the electronic spin becomes close to a good quantum 
number and Fig. 4(a) shows the energy levels split into 
two branches each with a different electron spin. The 
energy levels in each branch look like they move together 
as the applied field is changed. Figure 4(b) shows that 
this is only roughly the case and what look like straight 
parallel lines do have some curvature and structure. It is 
in this area where the weakly mixed ZEFOZ transitions 
occur. 

The transitions in the weakly mixed regime can be ex- 
pected to provide the longest coherence times as in gen- 



eral they are of lower curvature than those in the strongly 
mixed regime. The transition strength of the strongly 
mixed transitions will be much larger as the ZEFOZ tran- 
sitions in the weakly mixed regime are essentially nuclear 
spin transitions. Even though the strongly mixed regime 
is characterised by small Zeeman splittings compared to 
hyperfine splittings, a reasonable number of these tran- 
sitions occur at relatively high fields (>1 kG). This is 
because g is highly anisotropic, meaning that there are 
certain field directions for which /i#B • g • S is of com- 
parable magnitude to I • A • S resulting in energy levels 
which are still highly mixed. 

To obtain an estimate of the coherence times (T2) that 
should be attainable in Er 3+ :Y2SiOs the dephasing can 
be simply modelled as 

-L=S 2 (A£) 2 , (16) 

J-2 

where 62 is the curvature of the ZEFOZ transition and 
(AB) 2 is the variance of the magnetic field fluctuations. 
From reference [4] the magnetic field fluctuations expe- 
rienced by praseodymium ions in Y2SiOs are estimated 
as (AB) 2 = 0.0196 G 2 . For erbium dopants this could 
be considerably smaller due to the larger magnetic mo- 
ment of erbium inducing a much larger frozen core and 
so lifetime limited coherence times of the order of 200 ms 
are a possibility. In order to get the very long coher- 
ence times promised by the much flatter weakly mixed 
ZEFOZ transitions one needs to get around the problem 
of much faster spin lattice relaxation rates at high fields 
(26| . One way to solve this problem would be to operate 
at temperatures smaller than the splitting between the 
two branches and use a ZEFOZ transition in the lower 
branch. In this way the electron spin would always be in 
its lowest energy state and its relaxation rate would be 
immaterial. 



C. Transition Strengths and Hyperfine State 
Mixing 

More than 700 ZEFOZ transitions have been found in 
sites 1 and 2 of Er 3+ :Y2Si05; however, an important con- 
sideration is how easy it is to excite the transition either 
directly using an RF field or, of particular importance 
for quantum memories, via Raman transitions using the 
bottom level of the 4 ii3/2 excited multiplet. Precise val- 
ues for the Raman transition strengths are not currently 
available because the excited state spin Hamiltonian is 
not known. However, arguments about what transitions 
are likely to be strong and what are likely to be weak 
can still be made. For ZEFOZ transitions in the strongly 
mixed regime the mixing greatly relaxes the spin selection 
rules, allowing A-systems to be formed. This is especially 
the case around avoided crossings, which is where many 
of the ZEFOZ transitions occur. An illustration of this 
is provided by Baldit et al. where many A-systems were 
found at zero magnetic field (27| . 
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For the weakly mixed regime the electronic Zeeman 
term dominates and the rest can be treated as a pertur- 
bation. This means that ms is a good quantum number 
and the electron spin states are simply angular momen- 
tum eigenstates with the quantization direction given by 
B • g. For each value of ms there are eight hyperfine lev- 
els as the hyperfine contribution (A) is large compared to 
the nuclear quadrupole (Q), at least up to very high fields 
where the nuclear Zeeman term becomes important. By 
projecting the hyperfine term onto one of the ms = ±1/2 
branches it becomes of the form ±1 • A • (l/2n), where 
the unit vector n is the quantization direction for the 
electron spin. So like the electron spin, the dominant 
effect in the nuclear spin splitting for each of the two 
branches has the same form as a Zeeman splitting. This 
leads to nuclear spin states which are also angular mo- 
mentum eigenstates, but with a quantization axis that is 
in general different to that for the electron spin. So in the 
weakly mixed regime the transitions that leave the elec- 
tron spin unchanged will be allowed and can be driven 
with a RF magnetic field so long as the usual Ami = ±1 
condition is satisfied. 

The nuclear spin contribution to the transition 
strength of an optical transition is simply the overlap 
between the ground state and excited state nuclear spin 
states. If the quantization axes are the same for the nu- 
clear spins in the ground and excited states then Raman 
transitions will not be possible because the nuclear spin 
states will only overlap if the two mj values are the same. 
However, because the electronic Zeeman (g) and hyper- 
fine (A) tensors are anisotropic, and are different for the 
ground and excited states, the quantization axes will be 
different. The extent to which a Raman transition can 
be found will depend on the angle between the quantiza- 
tion axes. For small angle differences the change of basis 
matrix between the ground state and excited state en- 
ergy levels will be clustered around the diagonal, making 
Raman transitions possible only for small values of Am/. 

Thus as an alternative to calculating the optical tran- 
sition strengths, which would require knowledge of the 
spin Hamiltonian for the excited state, the values for 
| Ami | are given. Whilst this does not explicitly deter- 
mine whether it will be possible to create a A-system 
from the transition, it does give some indication. If Ami 
is small the chances of finding a strong A-system are much 
higher. 

The strengths of the ZEFOZ transitions, when driven 
with an RF field, can be calculated directly from the spin 
Hamiltonian using the probability amplitude method 
[30| . It is assumed that the transition is being driven 
by an AC magnetic field (Bac) oscillating at uj (where 
uj is the transition frequency). This situation can be de- 
scribed by the interaction picture Hamiltonian: 

H = H +H h (17) 

where Hq is the spin Hamiltonian of Eq. (|5]) and can be 



written as: 

H = (\a) (a\ + \b) (b\)H (\a) (a\ + |6) (b\) 

= hcu a \a)(a\ + hcu b \b)(b\, (18) 

where uo ai uo b are the energies of the levels involved in the 
transition. H\ represents the interaction of the atoms 
with the field: 

Hi = B AC COs(ut)[fl B g • S + UnQrJ] 

= B AC cos(ujt)(\a} (a\ + \b) (b\)[fi B g • S 

+ fJL n g n I\(\a) (a\ + \b)(b\) 
= B AC cos(ujt)[p aa \a) (a\ + p ab \a) (b\ 

+ Pba\b)(a\+p bb \b)(b\}. (19) 

Pab — (a> \ PBg ' S + Pn9ri\ \ b) is the magnetic dipole mo- 
ment of the transition from \a) to \b). 

As a measure of the transition strength the magnetic 
dipole moment of the transitions has been calculated. 
The magnitude of this is listed in Table IIVI for spe- 
cific ZEFOZ transitions. The ZEFOZ transitions with 
the lowest curvatures tend to have weaker transition 
strengths than those with higher curvatures. However, 
by making a compromise between curvature and transi- 
tion strength, several transitions have been found that 
have the potential to enable measurement of extended 
coherence times. Table llVl lists the properties of ZEFOZ 
transitions with curvatures ~ 10 3 HzG -2 and transition 
strengths > 10 5 HzG -1 . 



V. CONCLUSION 

In this work the ZEFOZ transitions for optical tran- 
sitions in Pr 3+ :YAG and ground-state hyperfine transi- 
tions in Er 3+ :Y2Si05 were calculated. This shows that 
the ZEFOZ technique which has previously only been 
used on the hyperfine transitions of non-Kramers' ions is 
also applicable to optical transitions and also to Kramers' 
ions. 

In Pr 3+ :YAG the terms of the excited-state enhanced 
nuclear Zeeman tensor were measured using Raman- 
heterodyne spectroscopy. The second-order Zeeman ten- 
sor was then calculated for the ground and excited states. 
Based on these parameters and those measured by Wang, 
the location of the ZEFOZ transitions were determined. 

In Er 3+ :Y 2 Si0 5 the location of ZEFOZ transitions for 
both sites were calculated based on the spin Hamilto- 
nian parameters measured by Guillot-Noel et al. Transi- 
tions were found for site 2 that have very low curvatures; 
therefore, it should be possible to significantly lengthen 
the hyperfine coherence time. The strength of the hy- 
perfine transitions was calculated, and it was determined 
that ZEFOZ transitions exist that simultaneously have 
low curvatures and high transition strengths. 
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TABLE IV: ZEFOZ transitions in Er 3+ :Y 2 SiO s where a compromise has been made between transition strength and maximum 
curvature. For all of these transitions |Araj| = 1. The levels are labelled 1 - 16, from lowest to highest frequency as depicted 
in Fig. SJa). 



Site Level 1 


Level 2 


A/fa vi mn im 

IVXCIj-A.1111 Li 111 

Curvature (HzG -2 ) 


B Dl (kG) B D2 (kG) B fe (kG) 


|B| (kG) 


TVp n sit inn 

11 CLllOl u 1W11 

Strength (kHzG -1 ) 


TVp n sit inn 

11 CLllOlulWll 

Frequency (MHz) 


1 7 


8 


1193 


11.43 


-6.136 


13.73 


18.89 


279.1 


905.7 


2 10 


11 


225.7 


-7.992 


3.529 


43.36 


44.23 


361.9 


408.0 


2 7 


8 


392.7 


-13.30 


9.169 


88.78 


90.23 


243.0 


239.8 


2 11 


12 


1153 


-1.029 


-4.825 


0.2901 


4.942 


887.0 


324.0 


2 12 


13 


1458 


-2.264 


-6.259 


6.567 


9.350 


427.0 


105.9 


2 5 


6 


1615 


0.5843 


4.024 


1.173 


4.232 


1292 


322.2 


2 6 


7 


1790 


-0.5366 


0.2861 


2.748 


2.815 


2796 


183.8 


2 6 


7 


1805 


-0.2467 


-0.7551 


1.411 


1.619 


4956 


465.2 


2 12 


13 


1812 


-0.7285 


0.09411 


4.776 


4.832 


4276 


417.9 
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